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Single-site dynamical mean field theory is used to determine the magnetic and orbital-ordering 
phase diagram for a model of electrons moving on a lattice with three orbital states per site and with 
the fully rotationally invariant Slater-Kanamori on-site interactions. The model captures important 
aspects of the physics of transition metal oxides with partially filled t2g shells, and of electron-doped 
Ceo- We introduce an unbiased, computationally simple and inexpensive method for estimating the 
presence of two sublattice order, determine the regimes in which spatially uniform and two-sublattice 
spin and orbital orderings are present and give physical arguments for the origins of the different 
phases. Guidelines are determined for optimizing the presence of ferromagnetism, which may be 
desirable in applications. 
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I. INTRODUCTION 

Orbital degeneracy is believed to play a crucial role in 
the physics of many "strongly correlated" materialsi of 
current interest, including "early" transition metal oxide 
compounds such as the lanthanum/strontium titanatcs^ 
and vanadateS)^!^ "late" transition metal oxide com- 
pounds such as the Sr/Ca ruthenate serieSf^i^ and non- 
transition metal oxide compounds such as the ^nCgo 
seriesj^ An important aspect of the physics of these com- 
pounds is the interplay between electron itineracy and 
the rich multiplet structure arising from the projection 
of the Coulomb interaction onto the orbitally degenerate 
on-site subspace. 

The single-site dynamical mean field approximationiS 
has had great success in treating the physics of materi- 
als in which only a single orbital is relevant. However, 
technical complications associated with the proper treat- 
ment of the full multiplet interactions have until recently 
caused difficulties in the application of this formalism to 
the multiorbital case. The single-band applications were 
often based on the use of the Hirsch-Fye quantum Monte 
Carlo mcthodii which does not have a straightforward 
generalization to the interactions present in the multior- 
bital case. A variant of the Hirsch-Fye method involving 
multiple auxiliary fields has been developed; however this 
method encounters a severe sign problem and difficulties 
with preserving rotational invariance in practiceJ^ Most 
published studies therefore involved approximations in 
which the exchange and pair hopping terms of the full 
Coulomb interaction were either neglected, treated via 
uncontrolled analytical approximation a ^^d'^d^ or approx- 
imated in a way which breaks the rotational invarianceii^ 

The recent development of continuous-time quantum 
Monte Carlo method a^^d^d^ along with improvements 
in the exact diagonalization technique^'^'^ have made 
feasible a comprehensive theoretical treatment of mul- 
tiorbital models with the fully rotationally invariant 



interactions 

In this paper we use single site dynamical mean field 
techniques to determine the orbital and magnetic order- 
ing phase diagram of a three orbital model. The work 
is an extension of a previous study of the paramagnetic 
phases of the modeli^ We show that some but not all of 
the Mott transitions are preempted by magnetic and/or 
orbital ordering transitions, determine the relation of the 
magnetic phase diagram to the non-fermi-liquid instabil- 
ity discussed in previous work^^ and demonstrate that for 
the model we study ferromagnetism exists only at carrier 
concentrations n > 1 and for very strong correlations. 
On the technical side we introduce a computationally in- 
expensive method to estimate the location of the phase 
boundary separating phases with uniform or two sublat- 
tice orbital and magnetic order from phases with no long 
ranged order, and where the order occurs estimate the 
transition temperatures. 

The rest of this paper is organized as follows. Sec- 
tion ini presents the model to be studied and the meth- 
ods we use, including the simple approach for identifying 
phase boundaries. Section [Till presents the phase diagram 
and gives physical arguments elucidating the origin of the 
various ordered phases, Section Hvl presents some results 
on the temperature dependence of phase boundaries and 
Section |V] is a summary and conclusion, outlining the 
implications of our results for experiments and prospects 
for future work. 



II. MODEL AND METHODS 
A. Model 

We analyze the "three band" model defined by the 
Hamiltonian 



H — i^band + ^int 



(1) 
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with 

^band 



(2) 



(ij)aP(j 



f3aa' '^iaa'^iPa' 



+h.c.). (3) 

^^band includes the usual electronic hopping between or- 
bital a = 1,2,3 on site i and orbital /? on site j, and 
the chemical potential /x. The fourier transform of t,^^ 
is a dispersion E(p) which is a matrix in orbital space. 
The two terms in Eq. ([3]) are the on-site interaction. 
The form of this term follows from symmetry considera- 
tions. The various Uapaa' terms are equal to (i) U when 
a = /3, cr ^ ct', (ii) U' = U-2J when P,a^a' and 
(iii) U'-J = U-3J when a 7^ /3, cr = cr'; where U{U') ac- 
counts for intra(inter)-orbital Coulomb repulsion among 
orbitals a and /3. J in the last term is the coefficient of 
the spin exchange and pair hopping terms. We shall be 
interested in cases in which the point group symmetry 
of the lattice guarantees that (in the absence of spon- 
taneous symmetry breaking) the on-site Green function 
G"^(i? = 0,0;) ~ 6ai3', this condition is satisfied in the 
materials listed in the Introduction. 



The self energies He.o are obtained from the solution 
of quantum impurity models (one for each sublattice A) 
of the form 



H, 



QI;\ 



a.cj 



hyb 



ifb 



ath; 



(6) 

with na^a,x the density of electrons of spin a on orbital a 
in the model pertaining to sublattice A and fi and A en- 
coding the chemical potential and any ligand fields aris- 
ing from the explicit breaking of the point group sym- 
metry. Spontaneous symmetry breaking is signaled by 
the appearance of a spin, orbital or sublattice-dependent 
Hartree term in the self energy. The interaction term 
Hioc consists of the on-site interaction terms of the orig- 
inal model. The remaining terms, i?hyb and i?bath sue 
bilinear in fermion operators and produce a hybridiza- 
tion function F^(a;) which is a matrix in spin and orbital 
space and whose form is fixed by the self consistency con- 
dition that the Green function of the quantum impurity 
model pertaining to sublattice A, G^^'^, be equal to the 
A - A component of Gioc (Eq. ©): G^^'^ = Gfj^^. 

In this paper we employ the semicircular density of 
states p{e) = \J — corresponding to a Bethe 

lattice with infinite coordination number and a fully bi- 
partite hopping Hamiltonian (Egamc hi Eq. ([4]) = 0). For 
this case the self-consistency equation in the general two- 
sublattice case becomes 



B. Dynamical Mean Field Approximation 



(7) 



To solve Eq. ([T]) we use the single-site dynamical mean 
field methodic which neglects the momentum depen- 
dence of the self-energy. To study magnetic or orbital 
orderings which spontaneously break the translational 
symmetry of the lattice down to a lower symmetry char- 
acterized by several sublatticcs, one must in principle 
introduce a quantum impurity model for each sublattice 
and take the self energy for the lattice problem to be lo- 
cal, but different on each sublattice. For simplicity we 
focus here on two sublattice orderings. We may then as- 
sociate a sublattice index A = e or o to the self energy 
(which also may be a matrix in spin and orbital space), 
distinguish hopping between the same and different sub- 
lattices and write the lattice Green function as a two by 
two matrix in sublattice space 



G 1 = 



Esamo(p) - 5]e(tj) 



Ediff(p) 

Ediff(p) LO - Esamo(p) - So(tj) 

(4) 

We assume that in the absence of spontaneous symme- 
try breaking (i.e. if Se = So ~ ^"ftrg) ^^^^ dispersions 
Esamo,diff are such that the local Green function 

Gioc(^) = J{dp)G{p,iu) (5) 

is proportional to the unit matrix in orbital space. 



with J-'A,a.a the hybridization function for sublattice A 
with orbital a and spin a and GB,a,a- the local green 
function on the other sublattice.— The self consistency 
equation for the translation-symmetry unbroken case is 
obtained by setting Gb ^ Ga in Eq. ([7]) and not consid- 
ering the B sublattice. We solve the quantum impurity 
model using the methods introduced in Refs. [islfigl The 
self-consistency condition is solved by iteration. 

For comparison we have also solved the model in 
the Hartree Fock approximation using standard meth- 
ods. For this problem the Hartree-Fock approximation is 
equivalent to replacing the full dynamical mean field self 
energy by 



J2{{U- 2J) {np^) + {U- 3 J) (np,) } 



(8) 



with a the spin direction opposite to a and a, (3 label- 
ing orbitals. The expectation values are determined self- 
consistently. The solution of the Hartree-Fock equations 
is simplified by the observation that the system becomes 
fully spin polarized before orbital ordering occurs. 
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FIG. 1: Variation of the on-site orbital occupancy (n) and 
magnetization (m) with number of iterations of the dynami- 
cal mean field self consistency condition, computed for a semi- 
circular density of states at density n = 1 with l3t = 50, 
U/t = 16 and J = U/6. Panel (A): iteration number depen- 
dence of the magnetization (summed over all orbitals); satura- 
tion indicates a ferromagnetic state with moment m « 0.65n. 
Panel (B): iteration number dependence of occupancies of 
the three orbitals, obtained from a translation-symmetry- 
unbroken (Eq. (O with Gb — > Ga) dynamical mean field 
equation; oscillations of the occupancy indicate two sublat- 
tice orbital ordering. Panels (C) and (D): iteration number 
dependence of occupancies of the three orbitals for sublattice 
A (panel C) and sublattice B (panel D), obtained from the 
general self-consistent condition Eq. (O; saturation of occu- 
pancies to a two-sublattice orbital state is evident. Error bars 
have the same size as the data points. 



C. Method of determining two sublattice ordering 

Finding orbitally and magnetically ordered solutions 
involves solving several quantum impurity models; it is 
therefore computationally more expensive than studying 
the paramagnetic case; in addition one must (as in any 
mean field theory) choose which symmetry breaking to 



investigate. We have found, however, that the presence 
of two sublattice order is indicated, to a good approxima- 
tion and in an unbiased way, by an oscillatory behavior 
in the solution of the dynamical mean field equations 
using the self-consistency condition appropriate for the 
case with unbroken translational symmetry. In studies 
of the symmetry-unbroken phases of quantum impurity 
models it is common practice to enforce the lack of or- 
der and improve the statistical accuracy by symmetrizing 
the hybridization function at each iteration of the self- 
consistency equation. However, if this is not done, then 
the dynamical mean field self consistency procedure will 
sometimes fail to converge, exhibiting instead an oscil- 
lating behavior which, we argue, signals the presence of 
two sublattice ordering. 

Figure [1] presents a particular example: results of a 
dynamical mean field solution of the three orbital model 
with semicircular density of states, interaction strength 
U/t = 16, inverse temperature (3t = 50, and density 
n = 1. For these parameters the model is in its Mott 
insulating regime, and as will be seen has ferromagnetic 
spin ordering and two-sublattice orbital ordering. The 
dynamical mean field equations are solved by iteration 
from a non-ordered seed state and the results are plotted 
against iteration number. Figure [T]A presents the mag- 
netization summed over all orbitals (the result is essen- 
tially the same whether the full Eq. (O or its translation- 
symmetry- unbroken special case is used). One sees the 
magnetization increase and saturate. Figure [Tl3 presents 
results obtained for the occupancies of the individual or- 
bitals using the version of Eq. ([7]) in which the transla- 
tional symmetry breaking is not allowed (Gb Ga)- 
One sees that a stable solution is not obtained. Instead, 
there is an oscillation between two states with different 
orbital occupancies. To demonstrate the meaning of this 
oscillation we present in panels C and D results obtained 
from a solution of Eq. ([7]), which explicitly allows for 
two sublattice translational symmetry breaking. Panel 
C shows orbital occupancies corresponding to one sublat- 
tice and panel D to the other sublattice. Comparison to 
panel B shows that the two states between which the self- 
consistent equations oscillate correspond to the states of 
the two sublattices in the ordered solution. We have 
also investigated other parameters with similar results, 
although care is sometimes required because when the 
two sublattice equations are used with a symmetric initial 
condition, two sublattice order may sometimes require a 
large number of iterations to develop, either because of 
a small Lyapunov exponent describing the growth of the 
broken symmetry phase or because the symmetry unbro- 
ken state is locally stable, so that one must wait for a 
fluctuation which is large enough to activate the system 
over a barrier. 

Of course, the existence of a two-sublattice solution 
to the mean field equation which persists for many iter- 
ations indicates at best the presence of a locally stable 
phase, and the parameter at which such a solution ap- 
pears or disappears may indicate a spinodal point, as in 
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many circumstances the transitions between phases are 
first order. An energy analysis would be required to de- 
termine the precise location of the phase boundaries, and 
the energy difi^erences involved are often very small (for 
example, the energy of the Ol paramagnetic insulating 
state at n = 1 and U = 16t is found to differ from that 
of the paraorbital paramagnetic state by an amount of 
order 0.3% of the total energy). 

The phase boundaries presented in the rest of this pa- 
per are obtained by determining the onset of oscillations 
in the solution of the translation-symmetry-unbroken 
equations and should be understood with these caveats. 



III. PHASE DIAGRAM 
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FIG. 3: Internal energy E computed as a function of density 
per site n at I3t = 50, U/t = 16. Vertical lines give the 
positions of phase boundaries presented in Fig. [2] 



Figure [5] presents the phase diagram in the plane of 
Coulomb repulsion and band filling obtained using the 
methods described above with J fixed to be U/6 and 
at the low temperature f3t = 50 (temperature 1/200 of 
the bandwidth). The phase diagram exhibits insulating 
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FIG. 2: (Color online) Main panel: phase diagram of the 3 
band model with semicircular density of states at f3t = 50 and 
J = (7/6 in the plane of particle density n and interaction 
strength U. The vertical lines indicate the Mott insulating 
phases at integral values of n. The magnetic state is labeled 
by P (paramagnetic), F (ferromagnetic) and A (two sublat- 
tice antiferromagnetic) while the labels 0(N) denote the 3 
classes of orbital ordering discussed in the text. The heavy 
dashed line (orange on-line) gives the boundary of the non- 
fermi-liquid frozen-moment phase discovered in Ref.[25l. Inset: 
Hartree-Fock phase diagram for magnetic phases of the same 
model. Magnetic phase boundaries are indicated by solid lines 
and orbital ordering boundaries by dashed lines. 00 and 
OS stand for the orbitally-ordered and orbitally-symmetric 
phases respectively. All transitions are second order except 
the FM-AFM transition and the orbital ordering transitions 
at (7 > 12t and small n. 



phases at integer band fillings n = 1,2,3 and sufficiently 
strong interactions. A new feature of the present paper 
relative to earlier work is the explicit inclusion of orbital 
and magnetic ordering. Comparison to Fig. 2 of Ref. [l^ 
shows that the inclusion of ordering does not shift the 
position of the Mott lobes appreciably: the critical U at 
the tip of the n — 1,2 lobes is essentially unchanged. 

Away from integer doping the phases are metallic. The 
various magnetically and orbitally ordered phases and 
their phase boundaries are indicated, as is the boundary 
to a phase discovered in earlier work^^ that is charac- 
terized by frozen moments and a non-fermi-liquid self 
energy. We have made selective studies of other values 
of J finding little impact on the qualitative features of 
the phase diagram. For example, for n between 1 and 2, 
increasing J to C//4 shifts the magnetic phase boundary 
downwards in U by about O.bt. The inset presents the 
Hartree-Fock phase diagram, which is discussed below. 

The phases we find are apparently stable against phase 
separation. Figure [3] shows the energy computed as a 
function of density for U = 16t; the upward concavity 
required for stability against phase separation is evident. 

B. Magnetic phases 

In the density range < n < 1 we found no evidence 
for magnetic order in the range of interaction strength 
studied. Ferromagnetic phases do occur at n = 1 (al- 
though for the largest U studied the Curie temperature 
falls below (St = 50 so no order is detected in our calcu- 
lations). Ferromagnetic phases also occur for 1 < n < 2. 
At n = 2 any uniform or two-sublattice magnetic tran- 
sition is below the temperature /3t ~ 50 of this study. 
We suspect that at n = 2 the magnetic order is actually 
incommensurate or has a higher order commensurability 
as was found for intermediate couplings in other mod- 
els with non-integer number of electrons per orbital per 
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sitei^ Ferromagnetism is also found for a range of n be- 
tween 2 and 3 while a two-sublattice antiferromagnetic 
state exists for n = 3. 

Note that the perfect nesting of the model we study 
implies that at half filling {n ~ 3) the antiferromagnetic 
phase is the ground state for all U > 0. The small U 
limit of the antiferromagnetic phase marks the point at 
which the Neel temperature falls below the temperature 
T = t/5Q used in the computations. 

A frozen-moment phase was recently reported by two 
of us in the three-orbital modeli^ The boundary to this 
phase is shown as the dotted line (brown on-line) in 
Fig. [21 Our results indicate that the frozen-moment 
phase is preempted by the ferromagnetic phase for strong 
correlations but that the frozen moment phase (which 
presumably indicates a non-uniform magnetic state; ei- 
ther disordered or a long-period spiral) exists over a wide 
range of dopings at weaker correlations. 

Remarkably, in the DMFT calculations the ferro- 
magnetic phase appears to be confined to interaction 
strengths of the order of or larger than the critical value 
for Mott insulating behavior and to densities greater than 
or equal to 1. Changing the Hunds coupling J over the 
physically relevant range J < U/3 does not change the 
results significantly. The DMFT results are in sharp con- 
trast to the Hartree-Fock results shown in the inset of 
Fig. [2] which indicate a much wider range of ferromag- 
netism and orbital order, extending in particular to much 
lower U and to densities n < 1. Our results suggest that 
LDA+U band theory calculations, which are in essence 
a Hartree-Fock approximation to the strong on-site in- 
teractions, may severely overestimate the range in which 
ferromagnetism and orbital order occurs. 

The absence of ferromagnetism at carrier concentra- 
tions n < 1 is similar to results obtained for two- 
orbital models by Momoi and Kubo2^ using a dynamical 
mean field method with an exact diagonalization solver 
and Held and VoUhardtii^ using Hirsch-Fye QMC with 
only the Ising component of the Hunds interaction re- 
tained, and also to very recent Gutzwiller approximation 
calculations^ which however yield a much wider range 
of ferromagnetism at ?i > 1 than is found in our calcula- 
tions. 



C. Orbital Ordering 

We now turn to the complex orbital ordering phase di- 
agram revealed in Fig. [2l All the d-orbital orderings pre- 
dicted by our DMFT calculations are staggered, rather 
than homogeneous in space. The nature of the phases 
is explicated in Table HI but the information about oc- 
cupancies must be interpreted with care. For example, 
minimizing the interaction energy at density n = 1 leads 
to a two sublattice ordered phase. In this phase, one sub- 
lattice has one orbital (say orbital 1) occupied and the 
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TABLE I: Characterization of orbital orders. gives the den- 
sity of electrons of the a orbit in sublattice A, etc. 5 is the 
deviation of the total density n from the integer value 1 for 
Ol and 03a or 2 (03b). Arrow lengths indicate magnitudes 
of densities. 



other two (2 and 3) empty, while in the alternate sublat- 
tice orbital 1 is empty and orbitals 2 and 3 are half filled. 
We believe the correct physical interpretation of the half 
filled state is that it represents an incoherent superposi- 
tion of the state with orbital 2 filled and 3 empty and the 
reverse, so that in the lattice the state corresponds to a 
highly degenerate set of states in which half of the sites 
in one sublattice have orbital 2 occupied while the other 
half of the sites on this sublattice have orbital 3 occu- 
pied. Supporting evidence for this interpretation comes 
from the measurement of the (712713) equal-time correla- 
tion function. At carrier concentration n = 1 it is found 
to be very small (of order lO^"'). 

This behavior can also be understood by considering 
the strong coupling limit. If we assume a fully spin polar- 
ized ferromagnetic state then the strong coupling Hamil- 
tonian describing the orbital ordering is (to leading or- 
der in t/U) an antiferromagnetic three state Potts model 
with nearest-neighbor interactions. On the cubic lattice 
this model is known to have a low temperature phase 
with precisely this structure)^ It is remarkable that this 
nontrivial state is correctly identified by our simple "find 
the oscillations in the DMFT iteration" procedure. We 
observe that considering higher orders in t/U would lead 
to longer ranged interactions which would lift the degen- 
eracy, leading to ordered states characterized by larger 
unit cells beyond the scope of our calculation. We sug- 
gest that the four-sublattice state discussed for LaTiOs 
by Pavarini et al^ is of this type. 

Doping changes the nature of the phase, converting the 
second sublattice from gapped to metallic (albeit with a 
low fermi-liquid temperature) in which the low T phase 
is (on long timescales) a coherent combination with each 
of the two sites half-occupied. In this phase the equal 
time (712713) correlation function is larger, of order 10^^. 
This phase, which is the stable solution of the DMFT 
equations, has within our 0.1% accuracy the same energy 
as a fully orbitally symmetric weakly ferromagnetic metal 
phase and it is possible that this latter phase (which we 
typically find to be unstable to the 01 phase) is the true 
ground state for some dopings. 
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FIG. 4: Magnetization, m, plotted against interaction 
strength JJ for different inverse temperatures /3 at density 
n — \ indicating a very rapid onset of magnetism as the mate- 
rial enters the Mott insulating phase and a weakly interaction- 
strength dependent Curie temperature at large U . 
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FIG. 5: Square of the magnetization (normalized to density) 
for various magnetic phases at fit — 50 and U/t = 16. the 
curve shown for density n = 1 corresponds to an insulating 
phase, while all other phases are metallic. The paraorbital 
phases (PO), n -"^ 1.2, 1.5, 1.75 have higher Curie temperature 
[kBT/t ~ 0.05 — 0.07) than the orbitally ordering phases. 



As seen in Fig. [21 further doping produces a transition 
to the 03a state. We have verified that raising the tem- 
perature shifts the phase boundary between Ol and 03a 
near n = 1 to the right as expected if Ol is the high en- 
tropy Potts model state whereas 03a is a lower entropy 
state. Similar arguments also apply for the competing 
orders 02 and 03b. However, the situation near n = 2 
differs in two important respects from that near n = 1: 
the 02 phase seems not to survive even an infinitesimal 
doping, and the 02 phase is entirely paramagnetic. 

The trends in the phase boundaries suggest that at 
strong coupling the intermediate disordered phase sepa- 
rating the 03a and 036 phases will be eliminated. We 
observe that at 5 = 1/2 the 03a and 036 phases become 
identical, implying that for large enough U where the 
disordered phase vanishes, the 03a and 03b states are 
adiabatically connected. 

Finally, we observe that the Hartree-Fock approxima- 
tion (inset to Fig. [5]) predicts a transition to an orbitally 
ordered state. The Hartree-Fock transition to the or- 
bitally ordered phase occurs at relatively large U, such 
that (within Hartree Fock) the spins are fully polarized. 
The transition is found to be second order except at very 
large U > 12t and small n, where it becomes first or- 
der. Again we regard the qualitative difference between 
the Hartree Fock and DMFT results as indicating that 
caution is necessary in interpreting LDA+C/ calculations. 



IV. TEMPERATURE DEPENDENCE 

In this section we consider the temperature dependence 
of the order parameters and the nature of the thermally- 
driven phase transitions. Figure [2] indicates that in the 



Mott insulating phase at n = 1 and at our chosen tem- 
perature Pt = 50, the system exhibits a phase transition 
from the ferromagnetic to the paramagnetic phase as U 
is increased above U ^ 18t. To clarify the nature of this 
phase transition, we present in Fig. |4] the dependence 
of the magnetization on interaction strength for different 
temperatures. The rapid and apparently T-independent 
onset of the magnetization at U ~ 11.5t suggests a first 
order magnetization onset which seems to coincide with 
the onset of the Mott insulating phase. More detailed 
studies are required to clarify the behavior near U = lit 
precisely. For the larger [/-behavior, comparison of the 
curves corresponding to different temperatures reveals a 
weakly U dependent Curie temperature, of the order of 
0.022t at [/ ~ U.5t, 0.02 at U ^ 18M and decreasing 
slowly for larger U, consistent with the ferromagnetic su- 
perexchange ^ t^ /U expected in this orbitally degenerate 
situation. 

To investigate the effect of doping we plot in Fig. [5] 
the temperature dependence of the square of the mag- 
netization, normalized to the particle density for vari- 
ous dopings between 1 and 2 at a C/ somewhat larger 
than the Mott critical values. At low T the spin po- 
larization is complete (to = n) except perhaps in the 
n = 1 insulating phase. One sees that metallic phases 
have higher Curie temperatures than insulating phases, 
and paraorbital phases higher Curie temperatures than 
orbitally ordered phases. The maximal value of the tran- 
sition temperature is seen to occur about half-way be- 
tween the Mott lobes, i.e. at n « 1.5, where the model 
is maximally metallic. This finding is consistent with 
the argument of Momoi and Kubo^ that the physics 
of ferromagnetism in multiorbital models is related to 
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V. CONCLUSIONS 
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FIG. 6: Relative orbital ordering strength for various phases 
at pt = 50 and U/t — 16. "Relative orbital ordering strength" 
is the deviation of the occupancy of one orbital from the 
paraorbital value (n/3) normalized to the value at our low- 
est measurement temperature ksT = 0.02t; mathematically 
[niikBT > 0.02t) - f)/(n^(fcsr = 0.02t) - f). The tran- 
sitions from the orbital selective to the paraorbital phase are 
apparently second order. The transition temperatures are not 
sensitive to the metallic (M/I), ferromagnetic (P/F) or order- 
ing (01/02/03) nature of the low T phase. 



the physics of 'double exchange': the strong interaction 
and non-vanishing Hunds couphng puts each site into its 
maximal spin state and the ferromagnetic transition tem- 
perature is then determined by the energetics of carriers 
hopping in the locally spin-polarized background. 

Figure \5\ also indicates that except very near to the 
Mott insulating phases, the temperature-driven ferro- 
magnetic transition is very steep, indeed within our ac- 
curacy apparently discontinuous. It is interesting that 
apparently first order behavior is most pronounced for 
dopings where there is no orbital order, while if orbital 
order is present the transitions seem more continuous. 

Figure [S] shows the temperature dependence of the 
order parameter for orbital ordering for several differ- 
ent carrier concentrations. The relative orbital ordering 
strength shown in the figure compares the deviation of 
from n/3 (the paraorbital value) to the corresponding 
value at low temperature (fc^T = 0.02t). The staggered 
orbital ordering has its critical temperature at around 
kBT/t ^ 0.04 — 0.05 and neither the value of the transi- 
tion temperature nor the order of the transition depends 
on whether the system is metallic or insulating, paramag- 
netic or magnetic. In the region where both orbital and 
magnetic ordering occurs the orbital ordering transition 
temperature is higher than the magnetic one. 



We have used the single-site dynamical mean field 
approximation to determine the magnetic and orbital- 
ordering phase diagram of a three-band model with re- 
alistic rotationally invariant multiplet interactions and 
have introduced a computationally simple and unbiased 
method for identifying the presence of two-sublattice long 
ranged order from an oscillation with iteration number in 
the symmetry- unbroken DMFT equations. We find that 
the Mott insulating phases at band fillings n — \ and 
71 = 2 per site are unstable towards a two-sublattice or- 
bital ordering of the degenerate type associated with the 
three-state Potts model on a bipartite lattice^^ while as 
expected at ?i = 3 the Mott phase is pre-empted by Neel 
ordering. It is interesting that the nontrivial ground state 
is the one identified by the unbiased method of examin- 
ing oscillations in the DMFT iterations. Ferromagnetism 
occurs at n = 1 and for n between 1 and 2 and for carrier 
concentrations greater than n = 2 but not too close to 
n = 3. The onset temperatures for orbital and magnetic 
ordering are low, of the order of 1-2% of the full band 
width. 

Remarkably, at n = 1 and n = 2 the orbital ordering 
is essentially co-terminus with the Mott phase: orbital 
ordering seems not to exist for interactions less than the 
Mott critical value, and only a tiny region of the Mott 
phases [U within a few percent of Uc) appears not to 
have orbital ordering. A high precision study of the be- 
havior for U very close to C/c2, beyond the scope of this 
work, would be required to verify the detailed behavior 
for U « Uc2- It is of course also possible that a longer pe- 
riod ordering (not considered in our work) would extend 
also into the metallic phases away from the Mott region. 
These issues, as well as the possibility of phase separa- 
tion, are important topics for further investigations. 

Ferromagnetism is found only for large interaction 
strengths, greater than or of the order of the critical val- 
ues needed to drive a Mott transition at n = 1,2 and 
only for carrier concentrations greater than n = 1 per 
site but not too close to n = 3 per site. As also dis- 
cussed by Momoi and Kubo, the physics of ferromag- 
netism in the multiorbital models appears to be related 
to the physics of 'double-exchange' and orbital selectiv- 
ity: at n > 1 and strong correlations one orbital becomes 
occupied by one (spin-polarized) electron; this occupied 
orbital acts as a 'core spin' whose orientation controls the 
hopping of the remaining electrons, leading to ferromag- 
netism similar to that in the manganites. For n < 1 the 
core spin effect is absent, and we expect that the physics 
of ferromagnetism is similar to that discussed in the con- 
text of the one-orbital Hubbard model^S where densities 
of states peaked in the band center (as is the case for 
the semicircle) disfavor ferromagnetism and densities of 
states peaked at the band edges favor it. We have found 
that varying the value of J within a physically reasonable 
range does not change the phase boundaries appreciably. 
However, the considerations from the one-orbital Hub- 
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bard model suggest that even for n > 1 variations in 
the density of states, in particular shifting the maximum 
away from the center of the band towards a band edge 
may widen the range over which ferromagnetism exists. 

The development of 'oxide heterostructures— has 
opened new avenues for material design. In the context 
of the LaTiOa/SrTiOa and LaAlOs/SrTiOa hctcrostruc- 
tures theoretical^SiS and experimental^ reports of fer- 
romagnetism have appeared. The materials of interest in 
these latter studies, however, were titanate-based, corre- 
sponding to densities n < 1. The theoretical predictions 
of ferromagnetism and orbital ordering were based on the 
LDA+U approximation, which is in effect a Hartree-Fock 
approximation to the local correlations. Our dynamical 
mean field results cast doubt on these theoretical predic- 
tions, and suggest that the magnetic behavior observed 



in Ref. [U may be due to something other than an in- 
terface ferromagnetic phase. Our results suggest that 
attempts to obtain a ferromagnetic electron gas at an 
oxide interface should focus on electron-doping a mate- 
rial with ?i = 1 or hole-doping a material with n = 2 and 
on arranging the bands to have an appropriate density of 
states. 

Three orbital models arise in a number of other phys- 
ically important contexts, including doped Cgo where, 
however, recent results suggest that lattice structure ef- 
fects may be important.— Extending our results to mod- 
els with more realistic densities of states is a high priority 
for future research. 
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